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Abstract 

We present an algorithm for simulation of SU(2) lattice gauge theory un- 
der the MA gauge and first numerical results for the theory without abelian 
monopoles. The results support the idea that nonperturbative interaction arises 
between monopoles and residual abelian field and the other interactions are per- 
turbative. It is shown that the Gribov region for the theory with the MA gauge 
fixed is non-connected. 



1 Introduction 

In the last ten years there have been extensive studies of abelian monopole dynam- 
ics inspired by conjectured dual-superconductor confinement mechanism [||, 0. The 
maximaly abelian (MA) gauge projection |^ was found having abelian dominance and 
monopole condensation in SU{2) gauge theory. It suggests the existence of an 
effective U{1) theory of confinement. 

But questions about role of Gribov's copies |^ and lattice artifacts are not clear 
yet. Calculation of monopole condensate like and direct investigation of the model 
without monopoles are also interesting. For these reasons a simulation of SU (2) gauge 
theory under the MA gauge constraints is desirable. 

This is the purpose of this work to present an algorithm for simulation of SU{2) 
gauge theory under the MA projection and first results of studying topology of the 
field manifold. Four- link analysis of the effective action for MA gauge is done in 0. 
We will keep here full effective action for the case of SU{2) lattice gluo-dynamics. 

In the next section we briefiy discuss the Faddeev-Popov operator for the MA pro- 
jection of SU{2) lattice gauge theory and a partial solution of the gauge constraints. 
We also give there a short description of a hybrid Monte Carlo algorithm of the simu- 
lation 1^. The numerical results and concluding remarks are presented in the last two 
sections. 
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2 SU{2) under the MA gauge. 



The MA projection for a lattice gauge field configuration {t/} in SU (2) theory is defined 
by i 

max R{{U} , {g}) = max J2 Trla^UiasUj^] 

where Ui = giUig'j are gauge transformed link variables and gi, gj are arbitrary SU{2) 
matrixes of the gauge transformation at sites i and j = i + fi. Written as a stationarity 
condition it reads 



5R 

or equivalently 



9=1 







where gi is an SU{2) matrix with g^^i taken to be zero. X = (JiX\ + (J2X2 + a^X^ is a 
traceless 2x2 antihermitian matrix and X-^ = aiXi + (72X2. 

Partition function for the theory with gauge conditions (|l|) could be written in the 
form 



5{X^) det{F) [dfi{U)] exp{-Sw{U}) 



(2) 



where dii{U) is the invariant integration measure, F is the Faddeev-Popov operator for 
the MA gauge and Sw is the Wilson action. 

According to Gribov's idea stationarity eqs.(|l]) should be complemented with 
the stability condition that the Faddeev-Popov operator be positively defined F > 0. 

On the lattice F = fij = — is a square symmetric matrix with nonzero diagonal 
elements [i = j) and elements for sites i and j connected by a link. On surface (|l]) /, 
reads 



X3,i 

X3, 



Tr{a2Uija2Utj) Tr{a2UijaiUi. 



Tr{aiUija2Utj) Tr{aiUijaiU^j) 

where Xi is defined by ([l|) and for Uij we suppose Uij = 
2N X 2A^-matrix where is number of sites. 



for i = j 



(3) 



for i ^ j 
Ujl- So F is a sparse real 
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2.1 Partial solution of the MA gauge constraints. 



Having in mind a hybrid scheme |T0[ to guide a MC simulation of theory with parti- 
tion function (Q) we now define a set of independent, with respect to constraints (|l]), 
variables {q}. 

Let a field configuration Cm a = {U} be a solution of (|I]). We consider three 
link variables S = {Usi,Uij,Ujk} of the configuration Cm a such that links Ls = 
{{s,i), {j,k)} form a continual path (3/-path) on the lattice. 

Equations ([I|), which include link variables of set S, read 



{Ul+a^U!, + ^asU!;)"- = {Ut^crsUs, + U.jCX.Ut,)^ (4) 

where the terms independent on S are omitted. Substitution Ug^ = Usi * exp {—ia^q) 
and f/J^ = exp {ia^v) * Ujk solves the first and last of equations (^). Two others are 
equivalent to 4 algebraic equations with 4 independent variables, say ui^ij , U2,ij , u^^ij 
and V, and allow to find a single solution and v{q). So, for any 3/-set S we can 
locally define a new set of variables 

Qs = {q, Qi, Qi, Qi, Qi}s = {q, Q^}s (5) 

where Qf- is the right-hand side of the i-th equation (^. If the link variables apart from 
S are kept fixed, then the gauge- fixing conditions require Q-^ to be constant, whereas 
g is a free parameter. 

Choosing a set of 3/-paths {S} such that Si nSj = ^ for i ^ j we can locally define 
a set of independent variables {q}. It allows to make transition Cm a Cm a where 
the both configurations lay on the surface defined by the gauge constraints (|l|). 



2.2 Hybrid Monte Carlo method for SU(2) gauge theory un- 
der the MA gauge. 

By following the well known procedure we can eliminate the determinant in (|^) and 
write the partition function in the form 

Z = J 5{X^) [dfi{U)d{<!>)] exp{-Sw{U} - $*F-i<l>) (6) 

where $ is an additional complex scalar field. It looks very similar to a partition 
function of gauge theory with pseudefermions except the (5-functions and form of the 
interaction matrix for the scalar field. 
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For numerical simulation of the system a Markov process should be constructed 
with the fixed-point distribution defined by and the only requirements for the 
transition probability of the process are detailed balance and ergodicity. 

The algorithm suggested in this paper includes two steps. First, we choose a set 
{S} of Ngite nonintersecting 3/-paths on the lattice. The choice is made at random 
from a uniform distribution. Second, we generate a new field configuration by hybrid 
Monte Carlo method [|l^, so the detailed balance and the gauge-fixing conditions are 



satisfied. We also repeat the transition for different {S}i to gain ergodicity. 

To specify the second step of the algorithm, let us take a configuration Cma{U} 
satisfying the constrains (|^) and select a set {S} of nonintersecting 3/-paths on the 
lattice. Restricted to the link variables associated with {S}, the desired probability 
distribution reads 



where Qq = Q and Jacobian A is given by 



Cma 



Tr{a,U,,a3U^)Tr[a3{U,,a3U^){Uika3Utk)] 



After integration over Q the probability distribution takes the form 



[dqd^] exp{-Sw{U{q)} - - ^ ln( A{Q, U) )) = [dqd^ exp(-5, 



(7) 



where the sum in the exponent is taken over all 3/-paths of set {S}. 

We apply the hybrid Monte Carlo algorithm pX| with evolution in pseudotime 
defined by Seff to generate a new MA configuration Cma- 

It is easy to see that Markov process converging to in variables Q yields a 
correct simulation of system (^. Moreover, if step-size in the integration of the equa- 
tions of motion is taken small, a quasicontinuous evolution in Langevin time keeps the 



configuration within the Gribov horizon [11 . 



3 Numerical simulation. 

The algorithm has been implemented for simulation of SU{2) gauge theory on the 
4^ lattice. We started with a configuration very close to zero field because matrix F 
for the completely ordered configuration has zero eigenvalues. The usual relaxation 
algorithm has been employed for initial MA gauge fixing. 

For the described in the previous section hybrid Monte Carlo algorithm, we chose 
duration 6t = 0.0005 for Langevin time steps and 400 steps for every run. The choice 
was not optimized but yielded acceptance rate ~ 0.7. The generated configurations 
were separated by 12 Langevin runs with different (random) sets of 256 3/-paths. It 
approximately reproduces a real number of degrees of freedom of the system. 
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The first 500 configurations were skipped for tliermolisation at [3 values 4.0, 3.0, 
2.5, 2.3 and ~ 2000 configurations aX j3 = 2.25. Tfie next 100 configurations were used 
to calculate average plaquette. Every configuration was searched for abelian monopoles 
1 12 1 but no monopoles were found. 

The results are summarized in Table 1. Here we also presented data on a one- loop 
calculations of the plaquettes [0 for 8^ lattice and values of the plaquette calculated by 
the heatbath algorithm. From the collected data we could observe that the calculations 
under the MA gauge are in an excellent agreement with the one-loop perturbative 
results at all (3, whereas the heatbath data essentially differ at /3 < 2.5. 



Table 1. 

Average plaquette in MA without monopoles, 1-loop expantion and full SU{2). 





^mon 


< PIma > 


P^one—loop 


< PIhb > 


2.25 





0.642(2) 


0.637 


0.593 


2.3 





0.647(3) 


0.645 


0.609 


2.5 





0.672(2) 


0.676 


0.655 


3.0 





0.729(3) 


0.733 


0.724 


4.0 





0.799(2) 


0.803 


0.800 



We additionally tested the stability of the generated configurations with respect 
to a small random gauge transformation with a further gauge fixing by the relaxation 
algorithm. In all cases the configurations appear to be stable. The situation changes 
crucially if matrix F is replaced by F with fa = fa and fij = for z 7^ j. At /5 = 2.25 
after ~ 200 sweeps starting from an almost ordered MA configuration we receive a 
nonstable one sometimes having abelian monopoles. So we come to the conclusion 
that it is the Faddev-Popov determinant that keeps us within the Gribov horizon. 



4 Conclusion 

We have presented an algorithm for the numerical simulation of SU (2) gauge theory 
under the MA gauge. The first numerical results give evidence that Gribov region for 
the MA gauge projection is split into at least two separated sectors. The configurations 
laying within the locus of the first zeros of det{F) have no abelian monopoles. Without 
monopoles the smallest Wilson loop - plaquette becomes very close to its one-loop 
perturbative value. It supports the idea that an effective theory could be built with 
nonperturbative interaction of the residual abelian field with monopoles and the other 
interactions are perturbative. 

This work is supported in part by ISF Grant of Long Term Research Program (No. 
RFWOOO). 
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